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A new criterion for pressure-driven interchange instabilities in cylindrical geometry is derived, 
based on an alternate use of the Energy Principle. This criterion is inequivalent to Suydam's 
criterion and does not contain the magnetic shear. In fact, it is shown that Suydam's criterion 
relates to the instability of the slow magnetosonic branch, while the present criterion relates to the 
Alfvenic one, which is the most dangerous of the two. These findings imply that pressure-driven 
modes often exist even if Suydam's criterion is satisfied by a large margin. 
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MHD instabilities are usually divided into pressure- 
driven and current-driven, depending on whether the 
destabilizing agent results from the interplay between the 
gas pressure gradient and the field-line curvature, or from 
the component of the current parallel to the magnetic 
field. The stability in cylindrical geometry with respect 
to both types of processes is an old and well-trodden sub- 
ject. It was actively studied in the 60s and 70s, both in its 
own sake and as an idealization of more complex geome- 
tries for fusion devices. Important general results have 
been established relatively early on, and are now part of 
the standard lore on the subject [j; among these Suy- 
dam's criterion 0] for interchange pressure-driven modes 
in an arbitrary screw pinch is directly related to the 
present work. Later developments have mostly been con- 
cerned with the stability properties of particular config- 
urations such as tokamaks or reversed-ficld pinches. 

In addition to fusion devices, astrophysical jets con- 
stitute another class of physical systems for which the 
existence of purely MHD unstable modes in cylindrical 
geometry is extremely relevant. However, there are only 
a handful of papers on current-driven instabilities in jets 
(0 Q an d references therein), and only a couple of 
papers on pressure-driven instabilities 6, 7]. This sur- 
prisingly small level of activity reflects the fact that the 
Kelvin-Helmholtz instability has long been thought to 
be the most important and dangerous for jet survival, a 
conjecture somewhat belied by recent numerical results 



Suydam's criterion stresses that the magnetic shear 
is an important stabilizing factor in cylindrical MHD 
columns. Indeed, although this criterion is only a nec- 
essary condition for stability, it is well-known that the 
magnetic shear plays an important role in the ability of 
nuclear fusion devices to survive the catastrophic disrup- 
tion usually produced by MHD instabilities. However, it 
has also long been noticed that pressure-driven unstable 
modes are easily found in MHD equilibria, with growth- 
rates nearly independent of their wavenumber, indepen- 
dently of whether Suydam's criterion is satisfied or not 
(see, e.g., Ref. ||). 



One of the objectives of this Letter is to provide a 
theoretical explanation for this fact, and to show that 
the role of pressure-driven interchange modes has proba- 
bly been underestimated in fusion physics. This is done 
through the derivation of a new general criterion for the 
stability of pressure-driven modes, and by a discussion 
of the physical meaning of this new criterion as well as 
Suydam's. The criterion is obtained from an alternate 
use of the MHD Energy Principle [J OH EH 

For static equilibria, the MHD perturbation equations 
can be written as 



- Po u 2 £ = -V8P* + ST = F(£), 



(1) 



where 5T is the perturbation of the magnetic tension, and 
<5P* the total pressure perturbation, and £ the fluid ele- 
ment displacement from equilibrium (a Fourier transform 
in time has been performed). Specializing to cylindrical 
equilibria, it is convenient to project this perturbation 
equation in the three directions e r , en, e±_ (i.e., the radial 
direction, the direction parallel to the unperturbed mag- 
netic field, and the remaining one required to define an 
orthonormal frame); furthermore the displacement can 
be chosen as a Fourier mode: 



€ = £( r ) ex P [i(wt — m6 — k z z) 



(2) 



This choice reflects the fact that the equilibrium quan- 
tities (the magnetic field line helical structure, the mag- 
netic field components and the gas pressure) depend only 
on radius. The equations for the parallel and perpendic- 
ular displacement £|| an d £ i can be solved exactly as a 
function of J 2 and £ r 0,02111: 



£||,_L = £||,±(^ 2 ,£r), 



(3) 



leading to an exact second order equation for the radial 
displacement 0,E1E3 : 



ar 



+ C*r£ r = F r (£ r ). (4) 
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In this expression, A and C* depend in a complicated way 
on the equilibrium quantities and on the mode frequency 

U). 

Let us define 



and 



(5) 



(6) 



where, for any choice of £ r , the remaining displacement 
components £|| and £j_ are specified by Eqs. and ui 2 = 
u 2 [£, r ] is chosen to be the formal solution of the implicit 
equation 



SW r 



(7) 



It can be shown from the standard form of the Energy 
Principle that, if there exists a total displacement such 
that SW r < 0, then the system is unstable; furthermore, 
uj 2 , as given by Eq. Q, also satisfies to 2 = SW/5K, where 
8W and SK are the usual potential and kinetic energy of 
the standard Energy Principle |l3j . 

Obviously, this reexpression of the Energy Principle 
is useful only if Eq. (0 can be solved for ui 2 , which is 
in general a rather daunting task. Fortunately, an in- 
teresting and useful simplification of F r is suggested by 
the ballooning asymptotic ordering. In particular, let us 
look for trial displacements where m, k z satisfy k±L Q ^> 1 
while k\\L <C 1 (L Q is the characteristic scale of equilib- 
rium gradients). In other words, one considers displace- 
ments with large wavenumbers, except in the direction 
of the unperturbed magnetic field. Because of the mag- 
netic shear, this imposes that the displacement must be 
nonvanishing only in the close vicinity to the magnetic 
resonance associated to m, k, and implicitly defined by 
fen = 0. Furthermore, let us make the assumption that 
one can find a trial radial displacement £ r such that the 
frequency uj 2 which is the solution of Eq. Q obeys the 
ballooning ordering, i.e. 
ir 



(^v 2 A /Llci/L 2 ) » 



„2 1,2 „2j,2 



where to 2 



c 2 )k\ is the fast magne- 



equilibrium azimuthal, vertical, and total magnetic field, 
respectively. The width of the region where the radial 
displacement is nonvanishing can be quantified from the 
constraint that kuL a is small enough (see below) with 
the help of this expression. For large k±L a , this width is 
narrow enough so that one can neglect the variations with 
radius of all quantities (except feii) in the corresponding 
region. 

From Eq. JHJ (and the approximations just described), 
the quantity C* reduces to 



pro: 
.13] 



C* 



2S| 
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B 2 
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(l + P)r 3 B 2 



(9) 



In this expression, the system of units is chosen such 
that fi a = 1, and the plasma p parameter is defined as 
(3 = c 2 /v 2 A: i.e. as the square of the ratio of the adiabatic 
sound speed (for an adiabatic equation of state) to the 
Alfven speed; this definition differs from the standard 
one by a factor 7/2 (7 is the adiabatic index). 

With these simplifications, Eq. (0 now reads 



^ / Po\ir? = / A 



dr 



dr + LO 2 / p \£r\ 2 dr, (10) 



where, to leading order, A is given by 



r k\ 



(11) 



Solving for u> 2 yields: 



f r r °^(\tir\ 2 + l/k ± 2 \dt r /dr\ 2 )dr' 



(12) 



where ui 2 is defined in Eq. @ . The quantity S appearing 
in Eq. (|12fl measures the width of the region where the 
radial displacement is nonvanishing. It must be chosen 
such that the constraint 



v\k 



2 - 2 kl < 



(13) 



tosonic frequency. This assumption is easily checked a is satisfied, with the help of Eq. © (if the displacement 
posteriori. 

With these approximations, feii can be expanded to first 
order with respect to the distance to the magnetic reso- 



is non zero in regions where v 2 A k 2 ,c 2 k 2 > \to 2 \, it pro- 
duces positive contributions to the potential energy 5W r 
|l3|). It is apparent that any displacement satisfying this 



nance r„, to read 



constraint makes lu 2 < provided that 



BgB z (r-r \ d\n\h\ 

where h = rB z /Bg is the pitch length (its logarithmic 
derivative is the magnetic shear), and Bg, B z , B Q are the 



^=^(-2^ c + ^)<0. (14) 

In Eq. lO, K c = e r-[(e||.V)e[|] = -B 2 /(B 2 r) is an alge- 
braic measure of the field line inverse radius of curvature, 
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while k p = (dp /dr)/p is a measure of the pressure gra- 
dient scale (an adiabatic equation of state is assumed). 
This reexpression of iv 2 allows us to recognise the first 
term as the usual destabilizing term due to the gas pres- 
sure and the field line curvature, while the second (sta- 
bilizing) term arises from the plasma compression. As a 
matter of fact, this criterion generalizes a result due to 
Kadomtsev for the saussage (m = 0) mode in a Z-pinch 

The reader can easily convince himself that the con- 
dition required earlier on u> 2 is satisfied for reasonable 
choices of £ r if ui 2 ^ 0, thereby insuring the self- 
consistency of the result. In the whole procedure, and 
in contrast to common usage, neither the incompressibil- 
ity condition nor marginal stability has been assumed 
(actually, the two conditions are known to be tied to one 
another). It is precisely because of this feature that it 
was possible to derive the criterion. Note that the usual 
reasoning according to which the most dangerous modes 
are incompressible is valid only if incompressible modes 
with finite growth rates can be found, which is not the 
case for pressure-driven modes in MHD static equilibria, 
as just pointed out. 

Within the ballooning asymptotic ordering, there is 
another well-known limiting case of direct interest to the 
purpose of this Letter, that make Eq. Q analytically 
tractable: namely, to 2 > v 2 A / L 2 ,c 2 / L 2 > v A k 2 ,c 2 k 2 > 

\uj 2 \ pa 0. This limit allows us to recover Suydam crite- 
rion. Indeed, in this limit, the fast variation of k\\ with 
radius dominates over all others, and after some algebra, 
Eq. igj reduces to 



dr 



9 k 2 P' 

-^ir = 0, 



(15) 



where A; = (k 2 + m 2 /r 2 ) 1 / 2 . With the help of Eq. ©, this 
equation gives the behavior of £ as a function of radius 
as £ oc [(r — r )/r ] p ; p is real when 



vBi 



8P' > 0, 



(16) 



p is complex otherwise. In this last case, it is possible 
to show, with the help of the Energy Principle, that the 
fluid is unstable (see, e.g., Ref-llj for details). Therefore, 
Eq. Ijl6|l (Suydam's criterion [2j) constitutes a necessary 
condition for stability. 

This stresses the fact that both the new instability cri- 
terion Eq. I|14|) . and Suydam's necessary stability crite- 
rion Eq. (|16fl follow from Eq. Q), but in different limits, 
and are not mutually exclusive. A deeper insight into the 
relation between these two criteria is obtained by consid- 
ering the low magnetic shear limit. In this WKB 
approximation can be made on £ r with radial wavenum- 
ber k r satisfying \k±r\ 3> \k r r\ 3> 1, so that the term 



involving the derivative of £ r can be neglected in Eq. Q . 
Assuming lu 2 » \lj 2 \, v 2 A k 2 , c 2 k 2 ,, Eq. @ reduces to the 
following dispersion relation 0, El 



u> 



kfi(v 2 A + v 2 SM )+uj 2 
k 2 (kj - 2(3k p k c )v 2 sm v a 



(17) 



where v 2 SM ~ (3v A /(l + ft) is the slow magnetosonic 
square speed. 

A necessary and sufficient condition of instability is 



k 2 < 2(3k p k c , 



(18) 



which requires in particular 2f3n p K, c > 0, a very-well 
known constraint for the existence of pressure-driven in- 
stabilities. Note that Eq. (|18(1 generalizes a criterion de- 
rived by Kadomtsev for the m = 1 mode [fj 0, El ■ In 
the limit fcy — > 0, there is a small and a large root to the 
dispersion relation Eq. I|17|l . which, to leading order in 
kuL a are approximately given by 



2 2/3k p k c v a 2 

' J SM 5 K \\ ' 



and 



(19) 



lo 2 + ~ to 2 , (20) 

Furthermore, in the limit n p , k c — > 0, the large and small 
roots reduce to the standard Alfven and slow magne- 
tosonic root, respectively. Only one of these roots is un- 
stable at a time, depending on the sign of ui 2 . 

Note that the large root w\ satisfies the ordering as- 
sumed in the derivation of the new criterion; furthermore, 
for displacements £ r satisfying the WKB approximation 
assumed in the derivation of Eq. I|17[). Eq. (|12[) gives back 
Eq. H2U|) . Reversely, the large root mode cannot satisfy 
the ordering needed to rederive Suydam's criterion, in 
contrast to small root modes. 

These arguments show that Suydam's criterion relates 
to the stability of the slow magnetosonic branch, while 
the newly derived criterion relates to the stability of the 
Alfvenic one, a conclusion that can also be reached by an 
analysis of the global modes behavior El- Of the two 
branches, the Alfvenic one the most dangerous, because 
it has the fastest growth rate, so that the new criterion 
is probably the more relevant of the two to ascertain the 
stability of MHD cylindrical equilibria. 

The practical usefulness of this new criterion also fol- 
lows from the fact that oj 2 < uj 2 /2 always provides a 
rather precise order of magnitude estimate of the unsta- 
ble pressure modes growth rates, even for low wavenum- 
bers (down to m = 1, in fact |l5|); this feature also ex- 
plains why the most unstable pressure-driven modes have 
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a growth rate nearly independent of their wavenumber. 
Quite surprisingly, the new criterion does not contain 
the magnetic shear. In fact, the magnetic shear only role 
is to limit the mode width through Eq. (fP3|) 13j; it is 
because of this feature, and not because of its influence 
on the growth rate, that the magnetic shear can limit the 
disruptive power of Alfvenic pressure-driven instabilities. 

Under which conditions can one expect Eq. (|14(l to be 
satisfied ? 

It is well-known that, in a Z-pinch equilibrium, ui\ > 
is a constraint which is difficult to realize in practice. In 
more realistic equilibria, however, one expects the longi- 
tudinal component of the equilibrium magnetic field to 
decrease with increasing radius, so that the gas pressure- 
gradient is relatively weaker in such a configuration than 
in a Z-pinch equilibrium, a feature which stabilizes the 
plasma, unless the pressure gradient occurs only on a 
rather small radial scale. In fact, this is precisely what 
happens in some of the most commonly considered toka- 
maks or RFP configurations. Furthermore, fusion de- 
vices are collisionless over the time-scales under which 
they usually operate; the simplest (albeit crude) way to 
account for this is to set the adiabatic index 7 = in 
the problem, which suppresses the compression stabiliz- 
ing term. Therefore, either for one or the other rea- 
son, pressure-driven interchange modes are necessarily 
present in fusion devices' plasmas, down to low azimuthal 
wave-numbers m > 1 (in any case, large wavenumber 
modes should be stabilized by finite Larmor radius ef- 
fects); furthermore, the growth rates are comparable to 
the Alfvcn time-scale once the plasma (3 is of the or- 
der of a significant fraction of unity. These features 
strongly suggest that the role of pressure-driven inter- 
change modes in the generation of the sawtooth oscil- 
lations and edge localized modes in tokamaks has been 
underestimated, a possibility which will be more fully ex- 
plored elsewhere. These modes might also play a role in 
triggering the edge turbulent transport in fusion devices; 
on the other hand, it is now well-known that a velocity 
shear can both stabilize pressure-driven modes and 
reduce the MHD-driven edge transport [l7j . A veloc- 
ity shear can also produce turbulent transport by itself 
through hydrodynamical processes, but the flow curva- 
ture and the Coriolis force can severely limit the induced 
turbulent transport [TsUli^. 

In what concerns astrophysical jets, most models stud- 
ied in the literature are cold, i.e., the equilibrium is not 
provided by the balance between the pressure gradient 
and the magnetic tension due to the azimuthal field, but 
by a balance between the inertial force due to rotation 



and the magnetic tension [2(1 l2ll |22(. Destabilization 
by the rotation inertial force is expected in this situation 
[3 , in spite of the stabilizing action of the Coriolis force, 
with features very similar to the ones described above. 
A more extensive study of rotation-driven instabilities 
in astrophysical jets will be presented elsewhere, while a 
comparison of usual pressure-driven instabilities with the 
Kelvin-Helmholtz instability is performed in Ref. [lfij |. 
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